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Microscopic theory of domain wall dynamics under electric current is reviewed. Domain 
wall is treated as rigid and planar. The spin-transfer torque and forces on the wall are derived 
based on the s-d exchange interaction between localized spins and conduction electrons, 
treating non-adiabaticity expressed by the gauge field perturbatively. Effect of spin relaxation 
is also studied. 
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1. Introduction 



The present information technology is based on electron transport and magnetism. Mag- 
netism has been most successful as high-density storages such as hard disks. In magnetic 
storages, read-out mechanism of the information had so far several significant developments. 
The oldest idea of detecting magnetic information would be to use Faraday's induction law by 
scanning a read head (a coil) on the stored magnetic bits. More efficient read-out mechanism 
was developed by using anisotropic magnetoresistance (AMR) effect. AMR is a resistivity 
dependent on the angle between the magnetization and the electric current, which arises from 
the coupling of magnetization and electrons' oribtal motion due to spin-orbit interaction. 1 The 
resistivity change is only of order of a few %, but AMR is more efficient than using Faraday's 
induction used in magnetic tape and hard disks in early days. Magnetic head with higher 
sensitivity was developed by use of GMR (giant magnetoresistance) effect in thin magnetic 
multilayers discovered in 1988. 2 ' 3 Strong magnetization dependence of the resistivity arises 
from the spin-dependent scattering of the electron at the interface between thin ferromag- 
netic layer and non-magnetic metallic layers. Quite recently GMR heads are being replaced 
by TMR (Tunneling MR) heads, where the non-magnetic layer is replaced by an insulating 
barrier. 4 These rapid developments of read-out mechanism by use of solid state systems made 
possible so far the rapid increase of recording density. 
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1.1 Current- driven domain wall motion 

In contrast, write-in mechanism in magnetic devices is still based on the knowledge of 19th 
century, Ampere's law. A novel mechanism of controlling magnetization by use of current but 
without referring to magnetic field was first considered by Berger in 1978. 5 He pointed out 
theoretically a possibility of driving a domain wall by current directly. Domain wall is a twisted 
structure of local spins as shown in Fig. I. 6,7 
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Fig. 1. (color online) Illustration of a Neel wall (left) and a Bloch wall (right). Also shown are the 
wall position X and the angle </>o between the wall magnetization and the easy plane, which are 
collective coordinates describing a rigid and one-dimensional (or planar) wall. A is the thickness 
of the wall. The axes shown are those in spin space (they do not necessarily coincide with spatial 
coordinate axes in the present systems without coupling between spin and real spaces). The current 
direction is perpendicular to the wall plane. 



The thickness of the wall, A, is an important parameter governing the coupling between 
the wall and conduction electrons. In most cases of 3d transition metals, A is 10 — lOOnm, and 
is much larger than the Fermi wavelength of the electron, hp -1 . The wall is therefore in the 
adiabatic limit, where the elctron spin can adiabatically follow the local spin direction as it 
passes through the wall. We consider in this paper a domain wall which is planar and rigid. 
In this case, the wall dynamics is described by two variables, its position X and 0o> average 
angle out of easy plane (x-z plane in both walls in Fig. I). 8,9 Our calculation applies to both 
Neel and Bloch walls, since there is no coupling between spin and real space in our model. 
The current direction is perpendicular to the wall plane. 

Berger considered a domain wall under electric current, and discussed that s-d type ex- 
change coupling between the local spin and conduction electron spin is dominant interaction 
that drives the wall under current in the case of thin film (e.g., thickness less than ~ O.l^m), 
where the effect of induced magnetic field can be neglected. 5 In 1984, 10,11 he studied the 
effect of the force arising from the reflection of conduction electron by domain wall caused 
by this exchange coupling. This force was associated with a wall mobility introduced phe- 
nomenologically. The effect was found to be small in most cases due to a very small reflection 
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probability because the wall thickness is usually large compared with Fermi wavelength. In 
1986, 12 he argued that the exchange interaction produces a torque which tends to cant the 
wall out of the easy plane (angle cpo). This torque was later found to push the wall by different 
mechanism from exchange force, and turned out to be dominant driving mechanism, 13 and is 
nowadays called spin transfer torque. Based on this idea, an experimental study was carried 
out in 1993 14 on a thin film of NigiFeig. There, domain wall velocity of 70 m/s was reported 
at the current density of 1.35 x 10 10 A/m 2 applied as a pulse of duration 0.14/is. Although the 
experiment was quite successful, experiments on better-controlled systems are now required. 

After these works, there was no significant development until 1996 in studies on current- 
driven domain wall motion, when Slonczewski 15 and Berger 16 independently developed a the- 
ory of magnetization switching by spin-transfer torque in thin film or pilar structures. This 
spin-transfer torque is essentially the same as the one Berger has discussed for domain wall. 13 
The pillar system considered there was intensively studied after the works by Slonczewski and 
Berger since such a system is expected to be applied to a memory devices like magnetoresis- 
tive ramdom access memory (MRAM) that opeartes without magnetic field. Current-induced 
domain wall motion is also expected to be useful as a possible MRAM, and intensive experi- 
mental studies have started. 

Recently experimental studies have been carried out on submicron-size wires and domain 
wall motion induced by current has been confirmed. 17-20 In all of the early experiments, the 
current density necessary for wall motion is high, of order of 10 11 — 10 12 A/m 2 . Measurement 
of domain wall velocity was carried out by Yamaguchi et al. 21 by observing wall displacement 
by use of MFM after each current pulse of strength 1.2 x 10 12 A/m 2 and duration of 5/is. The 
average velocity was found to be 2 ~ 6m/s. In those experimental works, there seemed to be 
a certain threshold value for domain wall motion, around 10 12 A/m 2 , and the average wall 
velocity were rather slow, less than lOm/s. 

1.2 Recent theories 

Those experiments motivated theoretical studies to look into the problem in more detail. 
Microscopic derivation of equation of motion of domain wall under current was carried out by 
Tatara and Kohno. 9 They considered a planar (one-dimensional) wall and described the wall by 
the two collective coordinates, X and 0o, i.e. within Slonczewski's description. 8 The variable X 
represents the position of the wall, and <po describes a tilt of the wall plane. Considering a small 
hard-axis anisotropy case, other deformation modes than 4>q (such as change of wall width) 
were neglected (rigid wall approximation). The equation of motion with respect to X and 4>q 
was derived including the effect of conduction electrons via the s-d exchange interaction. The 
electron carrying a current was treated by use of non-equillibrium (Keldysh) Green function. 
The equation of motion derived was found to be essentially the same as that obtained by 
Berger long ago, 10 ' 13 indicating his deep physical insights, but the effects of the spin-transfer 
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torque and reflection force (momentum transfer) were obtained without phenomenological 
assumptions and ambiguities for the first time. Based on the obtained equation of motion, 
the wall motion under steady current was studied. It was found that in the adiabatic limit, 
where the reflection force can be neglected, and if in the absence of spin relaxation, there is a 
threshold current determined by the hard-axis magnetic anisotropy energy, K± . Thus the wall 
is intrinsically pinned by the internal degree of freedom, (j>Q. At large current, however, the 
wall gets depinned and its velocity becomes proportional to spin current (spin polarization of 
the current flow), j s , as is required from the angular momentum conservation. 

Numerical simulation was performed based on an equation of motion of each local spin by 
including the spin-transfer torque term in the adiabatic limit. 22 The result was similar to the 
analytical (collective-coordinate) study, indicating the existence of threshold current. Motion 
of domain wall under magnetic field and spin-transfer troque was solved in Ref. 23 Later Zhang 
and Li 24 and Thiaville et al. 25 proposed to add a new torque term in the equation, which is 
perpendicular to the spin-transfer torque. After Thiaville et al., 25 we call this torque term as 
the [3 term. Zhang and Li argued that the (3 term arises from spin relaxation of conduction 
electrons. 24 Thus the phenomenological equation of motion of local spin under current reads 

ry . ri3 ft 

S = B cS xS + -SxS- —(j s • V)5 - ^-[S x (j s ■ V)S] + r na . (1) 

Here B c q is the effective field arising from spin Hamiltonian, and a represents damping. The 
equation S = B c q x S + |Sx S has been well-known as Landau-Lifshitz-Gilbert equation 
describing magntization dynamics in a magnetic field B c r . Effects of current are represented 
by other three terms in Eq. (1). The third term on the right hand side, (j s ■ V)<S, represents 
spin transfer torque, the fourth one is (3 term, and the last term r na denotes non-adiabatic 
term, 26 which has not been taken account in numerical simulations. 

The [3 term turned out to modify the threshold current and the terminal velocity of the wall 
significantly when (3/a> i. 24 > 25 > 27 The threshold current in this case is determined by extrinsic 
pinning and the terminal wall velocity is determined by /3/a. 24 ' 25,27 Microscopic derivation of 
the /3-term was carried out by several authors. 28-31 Tserkovnyak et al. 28 calculated f3 based 
on a one-band model considering spin-relaxation of conduction electrons semi-classically and 
assuming spin dynamics of small amplitude. They considered the limit of weak ferromagneticm 
and found that f3 s { = a s {, namely, (3 due to spin flip is equal to the damping parameter caused 
by spin flip. They also mentioned that in general /3 s f 7^ a s f considering the effects of multiband 
or deviation from weak ferromagnetism. Their approach is, however, still phenomenological, 
treating the spin-flip process by a phenomenological spin-relaxation time in the equation of 
motion of spin. (3 S [ = a s f was suggested also by different phenomenological argument 32,33 (see 
also Ref. 34 ). Fully microscopic calculation of (3 and a due to spin relaxation was carried out 
on s-d model by Kohno et al. 29,30 using standard diagrammatic perturbation theory, where 
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effect of spin relaxation are taken into account consistently and fully quantum mechanically. 
The result indicated B { ^ a s {. The same result was obtained later in the functional Keldysh 
formalism by Duine et al. 31 Determination of (3 and a values needs a careful microscopic 
calculation, since they are quantities smaller by a factor of l/(e^r) compared with conventional 
transport coefficients. Phenomenological thermodynamic argument predicted (3 = a, 32 ' 33 but 
microscopic studies 29-31 indicate that it is wrong. The error would be because the argument 
of refs. 32,33 lacks consistent consideration of the work done by the electric current. 34 

Waintal and Viret 35 and Xiao, Zangwill and Stiles 36 studied the spatial distribution of the 
current-induced torque around a domain wall by solving Schrodinger equation and found a 
non-local oscillatory torque (r na in eq. (1)). This torque is due to the non-adiabaticity arising 
from the finite domain wall width, or in other words, from the fast- varying component of spin 
texture. The oscillation period is ~ Uf^ 1 (&f is the Fermi wavelength) and is of quantum origin 
similar to the RKKY oscillation. Ohe and Kramer 37 studied a wall motion solving the torque 
due to the exchange interaction numerically, including non-adiabaticity. Non-local oscillating 
torque was numerically studied by taking account of strong spin-orbit interaction based on 
Kohn-Luttinger Hamiltonian (i.e., in magnetic semiconductors). 38 It was shown there that 
the oscillating torque is asymmetric around domain wall and that this feature results in high 
wall velocity. Non-adiabaticity was studied further in Refs. 26 ' 39,40 

In this paper, we review recent developments in the theory of current-driven domain wall 
motion. We consider the case of a rigid one-dimensional (planar) domain wall. This rather 
drastic assumption turns out to be more or less valid when compared with the numerical 
simulation and some of the experimental data available at present. Detailed comparison of 
experimental results with the present study will also shed light on the role of deformation of 
the wall, which is the future target. 

2. Model 

For simplicity, we take a localized picture for ferromagnetism and adopt the s-d model, 
where the dominant part of the ferromagnetic moment is carried by localized d-spins, S(x,t), 
and they are coupled to conduction electrons via the s-d exchange interaction. Essentially 
the same description would hold for itinerant ferromagnets, where the ferromagnetic order 
parameter plays the role of S(x, t) above. 

The Lagrangian of the system consists of that of electrons L e , that of localized spins Lg, 
and the s-d exchange interaction H s d between them; 

L = L c + L$ — H sd . (2) 

Each term will be explained in the following. Starting from this Lagrangian, we will derive 
the equation of motion of a domain wall. Since the domain wall is a macroscopic object, we 
treat it classically, whereas conduction elctrons are treated quantum mechanically. 
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The electrons we consider are interacting with impurities (both non-magnetic and mag- 
netic) and with external eletric field. We denote the electron annihilation and creation oper- 
ators as c a (x,t) and ct{x,t), respectively, where a = db represents the spin state. The total 
electron Lagrangian is given by 

L e = ^2 c k (*7$f - £fc) Cfc - H imp - H s{ - H em , (3) 
k 

where e& = 4^ — €f and Cfc = (cfc+jCfc-) 1 (t denotes transposition). The spin- independent 
impurity scattering is described by 

Hj mp = VjmpZ qRic k+q Ck i (^) 

i fc,q 

where Vi mp represents the potential due to impurity and R{ represents position of random 
impurities. We approximate the potential as on-site, i.e., v- imp has no q-dependence. Treating 
the impurity scattering by Born approximation, H{ mp gives rise to a lifetime of the electron, 
r CT , whose inverse is given by 

= 2irn imp Vi mp i' (7 , (5) 

No- 
where u a is the density of states at the Fermi level and ni mp is the concentration of impurities. 

The density of states and lifetime are spin-dependent in general, but we neglect this spin 

dependence when we discuss the spin transfer and momentum transfer, to avoid unnecessary 

complexity. This spin-dependence becomes important and so will be retained later when the 

effect of spin relaxation is studied. 

The term H s { represents spin-flip scattering due to random magnetic impurities: 



flsf 



J d 3 x S'idix - H^ac)^ (6) 



where S 1 i represents impurity spin at site R[. A quenched average is taken for the impurity 
spin direction as well as for the impurity position. 

The interaction with electric field is expressed by use of charge current j and electromag- 
netic gauge field A em as 

H cm = - / d 3 xA era -j, (7) 



The gauge field is given by use of E = — A cm as 

A cm = i£-e iUot , (8) 

where E is the applied electric field which is spatially homogeneous, and f^o is its frequency 
to be set as J7q - ► at the last stage of the calculation. The total current is given in the 
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presence of gauge field by (e(< 0) is electron charge) 

J = — ^(hk - eA cm )c{c k . (9) 




Fig. 2. Polar coordinates <p) parametrizing the local spin direction, and unit vectors n, eg and e^. 



2. 2 Spin part 

We consider the magnetization, or local spins, of fixed magnitude S and varying direction 
n, and parametrize it by the polar coordinates (6,(p) as (Fig. 2.1) 

S = S(sin 9 cos 4>, sin 6 sin eft, cos 6) = Sn. (10) 

Deferring the effect of damping (friction) to the next subsection, the spin part of the La- 
grangian is given by 

L S = hsJ ^<XcOS 6 - 1) - flg, (11) 

where H$ is the Hamiltonian of local spin, which we will specify later. The first term is known 
as the 'kinetic potential', and describes the spin dynamics governed by a torque equation. It 
has the same form as the spin Berry phase in quantum mechanics, but here we treat localized 
spins as classical objects. In fact, the equation of motion is derived from Lg as 

S = 7 S S x S, (12) 

where -yB s = ^|r- is the effective magnetic field acting on localized spin (in the absence of 
conduction electrons). 

The meaning of the 'spin Berry phase' term can be understood if one note that the 
canonical structure is contained in this kinematical term in the Lagrangian. Let us demonstrate 
this within classical mechanics. The canonical momentum conjugate to 4> is defined as 

p <t> = —^- = ns z - hs. (13) 

5(f) 

Defining the Poisson bracket (times h) by {A, B} PB = {8A/d(j))(dB /8S z )-(dB /dc/>)(dA/dS z ), 
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we have {(f), S z }-p& = 1. By using S x i 2iS*y — \J — e~^~^ , we can derive the correct SU(2) 
algebra of the spin angular momentum as 

{Si, Sj}pB = £ijkSk- (14) 

The Hamiltonian of localized spin we consider is a general one with two anisotropy energies. 
The easy axis and a hard axis, chosen as z and y direction, respectively. We treat local spins 
in the continuum. The Hamiltonian is given by 

+ V pin , (15) 

where V p i a represents sample inhomogeneity leading to the pinning of a domain wall. For 
a wire of soft ferromagnet, the easy axis is in the wire direction to avoid surface magnetic 
charges (shape anisotropy), and a domain wall appears as the Neel wall (Fig. 1.1). The case of 
Bloch wall such as realized in a film with perpendicular magneic anisotropy is also described 
by H$. As for the spin-transfer and momentum-transfer processes, both types of domain walls 
show the same dynamics if the spin-orbit coupling is neglected in the electron system. 

2. 3 Damping 

In spin dynamics, damping (friction) plays an essential role. We know that the magneti- 
zation will eventually point to the direction of the effective field, B s . Simple torque equation, 
S = jB s x S (7 = ^ is gyromagnetic ratio), or h = ^B s x n, however, predicts only a 
precession around B s . This point was remedied by Landau and Lifshitz (LL) by adding a 
perpendicular torque, 

n = jB s x n + a LL (n x (n x B s )), (16) 

where the last term describes a damping torque, which tends to align n along B s . Gilbert 
later proposed another form of damping, which contains n, 

n = jB s x n — ao(n x h), (17) 

and this equation (17) is called Landau-Lifshitz-Gilbert (LLG) equation. (In the above, qll 
and ao are dimensionless parameters.) These two equations are essentially equivalent, and 
they describe correctly the decay of precession and the relaxation to the equillibrium direction, 
n || B. As we will see, damping terms of Gilbert form can be derived by integrating out the 
environment. (But damping torque has also higher order derivatives, and so both LL and LLG 
equations are approximations in the linear order of time derivative.) 

The damping term cannot be introduced in the Lagrangian without environment, as is al- 
ways the case for dissipation processes. We here treat the Gilbert damping by use of Rayleigh's 
method in classical mechanics, 42 by considering a quantity describing energy dissipation, 

W S ee ^-hS J ^r(0 2 + sin 2 ej> 2 ). (18) 



J 



(VS) 2 



K 



(ft 



) 2 + ^(Sy) 2 
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The equation of motion with damping included is given by 

d 5L S 5L S 5W S 

dt 5q 5q 5q 
where q represents 6 and (p. The derived equation is the LLG equation, (17). 



(19) 



2-4 Exchange interaction 

The most important interaction in our problem is the s-d exchange interaction 

H sd = -J sd j d 3 x S ■ (cW) = —M j d 3 xn- (cW), (20) 

between local spin and conduction electrons. Here M = J s dS is half the exchange splitting. 
An important point is that J sd is rather strong in 3d ferromagnets: J s d/ e F ^ 0(0.1 — 1). 
These values are indicated from experimental observations of large magnetoresistances such 
as GMR. 

Most non-trivial part of the theory is the treatment of this strong exchange interaction 
when local spin has a spatial structure and/or is dynamical. Fortunately, spin structures 
in 3d ferromagnets are slowly varying compared to the scale of conduction electrons. This 
is a conseqence of strong exchange interaction, J, between local spins, which is of order of 
1000K as indicated by high critical temperature of 3d ferromagnets (For Fe, Tc ~ 1043K). 
(Correctly, the typical length scale, A, is determined by the ratio of exchange energy and 
magnetic anisotropy (Eq. (32)).) Since many local spins within the scale of A are coupled, the 
spin structure is (semi-) macroscopic and its time scale is slow compared to that of electrons. 
From these considerations, the electron can go through the spin structures adiabatically. The 
condition for the adiabaticity can be given by a few different small parameters. The first one, 
introduced by Stern 43 in disordered case, 

H/{Mt) « 1, (21) 

justifies the perturbative treatment of non-adiabaticity by using a gauge field. The second 

small parameter expressing spatially slow variation is l/(k F X) <C 1. For spin transport in the 

absence of disorder, this condition would be modified to be 

1 ep 1 

If ~ <r 1 (22) 

k F X M ~ (k F+ - k F _)X ' 1 ' 

where the left hand side is a ratio of the precession time of conduction electron due to the 
exchange interaction, h/M, to the time needed for the electron to pass through the spin 
structure, X/v F , as proposed by Waintal and Viret. 35 

3. Gauge Transformation 

Under the condition of Eq.(21), the electron spin is polarized almost along the local spin 
direction. By use of local gauge transformation in spin space, such electrons are mapped to 
electrons in a uniform ferromagnetic state interacting with a gauge field. The local gauge 
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Fig. 3. Gauge transformation of electrons in inhomogeneous field (due to local spins, S) into electrons 
in homogeneous field produces an SU(2) gauge field, A. 



transformation is to choose the electron spin quantization axis along S(x,t) at each point 
(Fig. 3). 41 The deviation from perfect adiabaticity is described by an SU(2) gauge field, which 
is small and we treat it perturbatively. A new electron operator a = (a+, a_) is defined as 



c(x,t) = U(x,t)a(x,t), 
where U is a 2 x 2 matrix which we take here as 

U(x, t) = m • <r, 

with a unit vector m given by 



m 



sm — cos 0, sin — sin 0, cos — 
2 2 2 



(23) 
(24) 
(25) 



The derivative of the original electron is written as d^c = U(dn + iA^a in terms of new 



electron, where A, 



-iC/(a;,t)- 1 M ?7(£c,t) = A a a a is the SU(2) gauge field with A? 



(m x d^m) a (summation over a = x,y,z is suppressed). The gauge field is related to the 
derivative n, as d^n = 2A fJi x n. 

The free-electron part of the Lagrangian is written in terms of the a-electron as 

^ ci a (iftd t - e k )c ka = ^2 a ka( ihd t ~ e k )a k(7 - H A (26) 

fc,CT k,cr 

where Ha describes the interaction with spatial and temporal variation of local spins, ex- 
pressed by A^; 



k,q 



(jJk + 1) ■ A«(-q)) a{ +q a a a k 



(27) 
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for fjL = x,y, z, t. 

The total electric current, eq.(9), is modified by the SU(2) gauge field as 



e 

m 



(hk - eA cm )a ] k a k + A a (q)a ] k+q a a a k 



q,a 



and the interaction with external electric field is given by 



ha k a k + A i(<l)al +q (T a a k 



q,a 



(28) 



(29) 



(30) 



up to 0(E). To summarize the electron part, the Lagrangian is now written as L c — H sc i = 
L° - H A - H cm - H s{ , where L° = Y.ka a ka( ihd t ~ e ka)o-ka ~ H imp: with e ka = e k - aM, 
defines free electrons under uniform magnetization and impurity potential. The interactions 
which we treat perturbatively are shown in Fig. 4. 




Fig. 4. Diagrammatic representaion of the coupling of conduction electrons to domain wall (Ha), 
to external electric field (H cm ), and to magnetic impurities (H s {). Solid lines are electron Green 
functions with selfenergy due to impurity scattering included, and wavy lines denote the SU(2) 
gauge field representing spin structure. 



4. Domain Wall 

We consider a planar (i.e., one-dimensional) domain wall as realized in a narrow wire, 
where the spin configuration changes only in the wire direction, which we choose as z direction 
of coordinate space. (This direction is y-direction of spin space in the Bloch wall case (Fig. 
1.1). Note that spin space and coordinate space do not necessarily coincide in the symmetric 
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system we consider.) The spin part of the Lagrangian L§ (without the pinning potential and 
the s-d exchange coupling) allows a static domain wall solution, 

z — X 

cos# = tanh — - — , (31) 
A 

with cj) = 0, and X is an arbitrary constant. Here we have introduced a length scale 

A - /J, (32) 

governing the spatial scale of magnetic structure in general, and gives the thickness of the 
domain wall. 

The domain wall considered above is called Neel wall (Fig. 1, left), where magentization 
is changing in the spatial z-direction, which coincides with the magnetic easy axis. Other 
type of wall, called Bloch wall (Fig. 1, right), is also possible, if the easy plane (zx-plane) is 
perpendicular to the wire direction, y. This difference of wall structure does not affect the 
electron transport nor the spin torque if the spin-orbit interaction is neglected. 



4-1 Collective coordinates of rigid ID wall 

To derive the equation of motion of a 'rigid' domain wall, we here consider the collective 
coordinate description. 44 This treatment and the results are essentially the same as the one 
considered by Slonczewski 6 ' 8 in the context of dynamics under magnetic field. 

The idea is to consider the constant X in Eq.(31) as a dynamical variable, X(t). Then the 
angle (p(z, t) can be excited, too, and so another collective variable 

Mt) = J ^sm 2 9 <p(z,t) (33) 

needs to be treated also as dynamical, 45 where cos 9q = tanh z or sin 9q = 

-l 

. This is because X and <po are canonically conjugate to each other, as in- 
dicated by the fact that the first term of Eq.(ll) takes the form ex Xfo. 

In the absence of sample inhomogeneity and driving force, X describes a gapless zero mode 
owing to the translational symmetry of the system. If the pinning potential Vq is present, the 
energy scale of X will be Vq. Similarly, the energy scale of the 0o-mode is given by K±. Since 
the energy gap of the spin- wave mode is ~ y/KK±, the modes described by X and cj>o are low 
energy compared to others if the following condition is satisfied; 

V < V KK ±> K ± « K. (34) 

In this case, the low-energy wall dynamics is described by the two variables, X and (j)Q. Oth- 
erwise, the pinning and/or K± leads to a deformation of the wall, whose description requires 
other variables than X and <pQ. The condition (34) gives a criterion that such deformations 
can be neglected. 

Precisely speaking, we need one more condition that there is no linear coupling of spin- 



cosh - "f^ 
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wave modes to X or cj)Q. In reality, when Vq and K± are finite, there arise such linear couplings, 
and the wall dynamics is not closed in X and 0o i n a strict sense. This is quite natural since 
the pinning and K± results in a deformation of the wall whose description requires other 
variables than X and (fro- However, the condition (34) also assures that such linear couplings 
are small. We assume the condition (34) in this paper. 

4-2 Domain wall Lagrangian 

From these considerations, the Lagrangian for low-energy dynamics of a rigid wall is given 
by using no = (9q, cj>o) in the Ls- The result is L = L w + L e , where 

L w = HNS [ ^-4> Q - ^j-S sin 2 O ) - Vp in [n ] 




+ M / d^xno-s. (35) 



Here s = c^crc is the electrons' spin-density operator, and N = 2\A/a 3 is the number of spins 
in the wall. (A is the crossectional area of the wire.) 

The equations of motion of the wall are now obtained simply by taking variations with 
respect to X and 4>o and taking the expectation value of s = (s). Using 

Ws = ^s(^ + ^) t (36) 



in eq. (19), they are obtained as 

X _X_ 

T ~ Ms 



+ a o— = T^( F + F Pin), (37) 



X - a X4>o = -jjr- s sin 2( Po + ^5 Te ^- ( 38 ) 



Here force and torque due to electrons are defined as 

5X 



F = -(^) = -M I d 3 xV z n -s, (39) 



d A x ( 5 -^L )xS = -M I d a x(n x s). (40) 



5S 

We note that this set of equations, (37) and (38), is essentially the same as those obtained 
by Berger. 10 ' 13 What is new and essential in the present theory is that we have formal but 
exact expressions of force and torque, which we can evaluate by a systematic diagrammatic 
method. 

Defining each component of s as (see Fig. 2 for the definition of eg and es) 

s = s e e e + + s z n, (41) 
force and torque are represented in terms of sg and as 

Fn = 2M J d 3 x (sgA* + s^j) (42) 
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r e = M J d 3 x (s^e e - sge^) , (43) 

where Af, = ■ A M , A^ = eg ■ A^. Clearly, the component s z parallel to the local spin S does 
not affect the dynamics of S. 

5. Calculation of Electron Spin Density 

Our task is now to evaluate the electron spin density s in the presence of spin structure 
(Vn) and current flow, or in the presence of spin dynamics (n). Since the electron state is 
better defined in gauge-transformed (rotated) frame, we define Green functions with respect 
to the rotated frame. We use non-equillibrium (or Keldysh) Green function defined on complex 
time plane, 46 

9ka(t, t') = -i (Tca ka {t)a{ a (t')^ , (44) 

where Tq represents path order on complex time plane (Fig. 5) and (• • •) denotes averaging 
over quantum states, random impurities and thermal averaging. This non-equillibrium Green 
function contains besides retarded and advanced Green functions the lesser (greater) Green 
function, 

SfarftO =*(aL(Oafc CT (t)), (45) 



which gives directly the informaiton on the particle (hole) number and is is most useful 
in calculating physical qunatities. The quantum state is defined with respect to L[?, with 



C 



3-^ 



Fig. 5. Keldysh contour in the complex time plane used in non-equillibrium Green functions. 



normal impurities taken account in the standard ladder approximation (i.e., neglecting small 
corrections of 0(1/(€ft))). Thus we have free retarded and lesser Green functions as 

= (46) 

where f{u) = ^pj- 

The exact Green function (detnoted by G), taking account of interaction, Ha, H em and 
H s f, satisfies the same Dyson equation as the conventional time-ordered and retarded Green 
functions (but with time defined on a complex plane). This Green function is still defined on 
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complex time, t,t'. To obtain physical quantities, we need to map i's onto real time. 46 

We consider a slow local spin dynamics and assume that the SU(2) gauge field has only 
zero frequency component, A°(q, Q) = 5n t oA^(q). This is justified when Or <C 1. We can 
easily evaluate the spin density in rotated frame, s, defined by 

s(x,t) = (afo-a) = -itr[a a G < (x,t,x,t)]. (48) 

The spin density in the original frame is given by 

s(x, t) = 2m(m ■ s) — s. (49) 




Fig. 6. Contributions to the electron spin density in the rotated frame, in zeroth- and linear-oder in 
the applied electric field E (denoted by and , respectively). 



6. Spin- Transfer Torque on Domain Wall 

We consider a rigid one-dimensional domain wall represented by (31) (33). The correpond- 
ing gauge field is given as (q is along z-direction) 



4(q) = 



u, 



(50) 



i 



is a form factor of the wall and A\ = 0. The electron spin polarization 



where U 1 = cosh(|A g ) 
around the domain wall is then obtained as 26 

7rA 



Sg(x) 



2 LAI 



iq(z-X)„ 



(51) 



s<f>(x) 



7rA 

2LM 



iq(z-X) 



Un 



eX 



(52) 
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where s = (k F+ — k F _)/(127r 2 ) is the equilibrium spin density of electrons, kp = (kp + +kp_)/2, 
and P = j s /j = 5n/n is the polarization of current with n = n + + ?i_ and 5n = n+ — n_. 
Dimensionless correlation functions are given by 

xi%) = ^£ P fj ^- r (53) 

sV ^ e k+q -e k ±2M 

xf{q) = £ - MS (e k+q -e k ± 2M) , (54) 

kp V k,± 2 

k,± 



fe,± 

x^ ± - 9 y% +9 - efc i2M). (56) 

As seen, an d (* = 0> -0 ar i se from the real and imaginary parts of the factor, 

r , respectively. 

Let us look into the adiabatic limit by putting q = in x's. The spin density then becomes 

sf\x) = ^sme a (z-X)sfo, (57) 

s * d, w=-2s siii9 °< 3 -- Y »(4 + ^)' < 58 > 

Results (51) and (52) indicate that non-adiabaticity (finite q contribution, xi ) basically 
exchanges roles of and (f>; non-adiabatic contribution from current induces a spin polarization 
in the ^-direction, and drives tilt of the wall. This is exactly what is expected in the presence 
of spin relaxation as argued by Zhang and Li, 24 so non-adiabaticity and spin relaxation have 
essentially the same effect on the spin structure. 

These features are cleary seen in the equation of motion. The torque on a wall is obtained 

as 



\ J d?>X j d3x ' sin -0 Q ( x ) 8in 0o (V) 

X ~(0)/ /\ is Is 

syXi {x - x ) - 2^ -^Xi '(x-x) 



Noting X\\l) = 1 + 0(q/kp) 2 , we obtain 



(59) 



'sX-^-) +o((kp\y 2 ). (60) 



2A V e 

Roles of electron spin polarization on wall dynamics is summarized in Fig. 7. 
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Fig. 7. (color online) Spin polarization of conduction electrons (denoted by small arrows) around a 
domain wall (large arrows), and its effect on the domain- wall dynamics. Top: Equilibrium polar- 
ization along the local spin does not affect the dynamics. Middle: Under current, the spin-transfer 
torque is induced by the polarization component out of the wall plane (in spin space), which 
induces rotation of local spins within the plane, hence the translational motion, X. Bottom: 4>q 
is induced by the component in the ^-direction in the wall plane. Such polarization arises from 
non-adiabatic 9 and spin-relaxation 24 processes. 



7. Spin Relaxation 

The presence of spin relaxation processes in the electron system produces a new type of 
torque 

a 3 S 

r s f = a si Sn x h + Atf-j^-fa x (J B ■ V)n), (61) 

in the adiabatic regime, as first recognized in ref. 24 The first term is the Gilbert damping, 
and the second term is a new type of current-induced torque which is orthogonal to the spin- 
transfer torque. (a s f and /3 s f are dimensionless coefficients.) In ref., 29 the authors adopted 
quenched magnetic impurities to simulate spin relaxation processes in a microscopic Hamil- 
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tonian, and obtained the result as 

Q sf = irn s u 2 



Full Paper 



'si 



nn s u 2 
M 



2S 2 z v + v- + S\ (v% + v 2 _ 



(S 2 ± + S 2 )D + + ^(S 2 ± -S 2 )^ 



(62) 
(63) 



in terms of the density of states, v±, or u± = u + ± z/_ (note that the notation is different 
from ref. 29 ), the concentration (n s ) and the scattering amplitudes (u s S±, u s S z ) of magnetic 
impurities, and the degree of spin polarization, P = j s /j, of the current. 



8. Force 



The concept of force may be generalized to arbitrary spin structures based on Eq. (42). 
There are several types of forces corresponding to each torque. 26 ' 47 In particular, we now have 
three kinds of current-induced forces: F = F refl + F ST + F@. The first one 

F IC& = eN cPs j, (64) 

is a (non-adiabatic) force due to electron reflection by spin structure, and is related to the 
resistivity 48, 49 

' ^|^(g)| 2 <5K CT )5( efc+? ,„ CT ), (65) 



k,q,cr 

due to the spin structure. The other ones, F and F@, are finite in the adiabatic limit, and 
are physically different from the reflection force, _F rcfl . They are calculated as 

F i ST =4e^/ d3x n • c^ n x d ^ ( 66 ) 

Ff = d 3 x a-V)n • Vin, (67) 

up to 0{{kp\)~ 2 ). For ferromagnetic films, the spin-transfer force F" ST has a topological 
meaning since the quantity, n v = ^- f d 2 xn ■ (d x n x d y n), defined in two dimensions is a 
topological number (integer or half integer depending on boundary conditions). The force 
F ST is, in fact, a back reaction of the Hall effect due to spin chirality, 50,51 and was derived 
by Thiele 52 and Berger, 12 then based on a general relation between force and the torque, 47 
and also in the context of magnetic vortex. 53 Note that the adiabatic force is included in the 
(adiabatic) spin-transfer torque, while the reflection force is not. The second term due to the 
/3-term is written as 24,25,47 

F P = ^T'/Ws, (68) 

for a one-dimensional spin texture, where 7' = f d 3 x(V z n) 2 . These forces are schematically 
illustrated in Fig. 8. 

There are also forces induced by spin dynamics (S) in the absence of current, F = F^ + 
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Fig. 8. (color online) Schematic illustration of three kinds of current-induced forces acting on a spin 
structure. Electron reflection pushes the structure along the current flow (F rcti ), and Hall effect 
due to spin chirality results in a force in the perpendicular direction (F ST ). Spin relaxation results 
in a force along the direction of current (F^). 



f(l) . 

F/ 0) = S -jfx n -(hxd i n), (69) 

F^=aaSj ^n-fyn (70) 

The first term comes from the "spin renormalization" torque (see the next section), and the 
second term comes from the Gilbert damping. 

9. Equation of Motion 

For a rigid one-dimensional wall, = — ^rr^o + o((^fA) -2 ) and F ST = 0. (For a vortex 
or a vortex wall, F ST is finite.) Combining Eqs. (37), (38), (60) and (64), we finally obtain 
the equation of motion under current as 

X a? 

ls4>o + a— = —(3j + /pi n , (71) 

3 

1S X - a\<f>o = v c sin 2^» + ^P J ~, (72) 

up to 0((/cfA) -2 ), where the damping now includes the contribution from spin relaxation (eq. 
(62)), a = a + a sf , 

= 4[§ e2nR ™ A + P(3sf ' (73) 

is the total force due to electric current including the effect of spin relaxation (eq. (63)), and 
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/pin = -i^X6{i - \X\), (75) 



has the dimension of velocity. Pinning force represented by 

2XV 

nstf 2 

is due to an (extrinsinc) pinning potential V p i n (X) = ^M W Q 2 (X 2 — £ 2 )#(£ — |^|), which we 
approximate as harmonic with range £, pinning frequency Q (M w = N/ (K±X 2 ) is the wall mass 
and Vq = |^0 2 ^ 2 ). The factor of 75 = (S + |) /S, arising from F (0) , indicates that electron 
spin polarization s/2 contributes to the magnetization of the wall (if close to the adiabatic 
limit). This natural result indicates consistency of our calculation. (In the equaitons of motion, 
we have neglected small non-adiabatic corrections of o((fcpA)~ 2 ). Correctly speaking, the 
reflection force is also the same order of o((fci?A) -2 ), 48 ' 49,54 but we retain this term since 
physical meaning is clear.) 

Most important part of our work is the derivation of this equation of motion. We see 
that it is essentially the same as the one argued by Berger 10 ' 13 (with additional terms in- 
troduced by Zhang and Li, 24 and Thiaville et al. 25 ), indicating his deep physical insights. 
The phenomenological arguments were quite useful in discussing the adiabatic limit, where 
angular momentum conservation (adiabatic spin-transfer torque) governs the dynamics. Once 
non-adiabaticity and spin relaxation come in, fully quantum mechanical calculation as we 
did is required. Of particular future interest would be the quantitative first-principle estima- 
tions of torques and forces including material parameters (such as spin-orbit interactions) and 
geometry (pinning) based on the present formulas (42) and (43). 

10. Solution 

Let us look into the solution of the equations of motion, which are given in terms of 
dimensionless parameters by 



d- t [X - a0 o J = sin20 o + Pj, (76) 

Of (0o + ctX) = -V X6(Z/\ -\X\)+ pj, (77) 

where t = tv c /X, X = X/X, Q = £lX/v c = ^r\J^, P '= = ^3 and V = \<h 2 = 

^ jj^ , and we approximated as 75 = 1. (For details of the solutions, see Ref. 27 ) 

10.1 Threshold current 

Behavior of threshold current depends on the extrinsic pinning. 
(I) Weak pinning regime Under small current, j < 1, 4>q remains small and the wall 
dynamics is well described by X only. This is defined as regime I. Linearizing the sine-term 
in Eq.(38) as sin 20 ~ 20, we eliminate to obtain 55 ' 56 

(l + a 2 )d?X + \d- t X + n 2 X = Fj, (78) 
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where 1/f = 2a (l + §F ) , and Fj = 2/3j is a dimensionless force due to current. We consider 
the case of steady current and weak damping; 2Qf > 1. A solution satisfying the initial 
condition, X(0) = and d^X(0) = Pj, is obtained as 

X(t) = ^ ( 1 - e~w ( cosSK + -^UsinO'i 



ft 2 V V 29,'f 

Pi i 

+ -J-e~w sinWt, (79) 

where O' = \JVL 2 — The threshold (depinning) current is determined by |X(j)|inax = £) 
which is given by 

*»> = £ 3 fi (A S (80) 

By balancing the pinning force and the force due to magnetic field, Vq is expressed by the 
depinning field B c as 

V = | gfi B B c ^. (81) 

(II) Intermediate regime This regime, j > 0(1), could be important for application since 
the threshold is not sensitive to sample irregularities. Depinning in this regime is described 
by cpQ. 9 The reason is that the effective mass of "0Q-particle" , given by l/Vo 45 (see Eq. (82)), 
becomes lighter than the corresponding mass of "X-particle" given by 1/K±, and so "<^q- 
particle" is a better variable to describe dynamics for strong pinning. By eliminating X from 
Eqs. (76) and (77), we obtain 

(1 + a 2 )dU + adrfo ( 2 cos 20 o + V 



+ V sin2<f) +jV P = 0. (82) 

Thus (3 does not affect the dynamics of <^o- (Correctly speaking, this feature is specific to the 
harmonic pinning potential, and anharmonicity introduces the dependence on (3.) In fact, the 
/3-term can be eliminated from the equations of motion if one rewrite Eq. (77) in terms of 
X' = X — |£j (i.e., it just shifts the stable point of X). Even in the case with anharmonicity, 
we have numerically checked that (3 does not lead to important modification in this regime. 

From Eq. (82), we see that the energy barrier for <po vanishes when j c ~ -P -1 , irrespective 
of the pinning strength. Once </>o escapes from the local minimum, its velocity is given from 
Eq. (82) as 9j0o — 77" -This corresponds, according to Eq. (77), to a maximum displacement, 
I max ~ of the wall. Unless the pinning is extremely strong, i.e., if olVq < 1, |X max | 
exceeds £/A, i.e., depinning of X occurs as soon as cpo is depinned. Thus the threshold is 
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Table I. Summary of threshold current for different strength of extrinsic pinning. 



Pinning 




Threshold 


I-a: Weak 


< O(l), /3 < O(O) 


j c oc VKl^o 


I-b: Weak 


n < o(i), p > o(fi) 


j'c oc Vo/ [3 


II: Intermediate 


o(i) < n < 0(0 


jc OC K_L 


III: Strong 


0(0 <o 


j c oc Vb/a 



roughly given by j c ~ P 1 , and is actually found numerically to be 

e AS 12 

i~ l,l W S - < 83 > 

This story was presented in ref., 9 but the estimate of threshold current there was j c = 1. 
The reason for this difference comes from (3. In the analysis of ref., 9 where (3 = was assumed, 
even if X escapes from the pinning potential for current j > 0.7 P^ 1 , the terminal velocity 
vanishes if j < P^ 1 owing to the intrinsic pinning effect (i.e., <po reaches a steady value and 
X becomes zero). On the other hand, if P 7^ 0, steady motion of X is possible as soon as X 
escapes from the pinning. This is the reason why the threshold curreent is different for (3 = 
and (3 ^ in the intermediate regime II (Fig. 9). 

(Ill) Strong pinning regime The above result X max ~ j^- indicates that for extremely 
strong pinning, Vq > a -1 , the wall is not always depinned even after <pQ escapes from a 
potential minimum due to K±. The depinning occurs at 

Jc ~ -J-y (84) 

as pointed out in ref. 9 

Results for the threshold current are summarized in table I. It is interesting that such a 
simple set of equation of motion results in so rich behaviors. 



10.2 Wall speed 

After depinning, the wall dynamics is describged by the equations of motion, (76) and 
(77), with Vq = 0. The solution can be obtained analytically (see Eq.(31) in Ref. 27 ). We see 
that the wall dynamics is quite different for j > j a and j < j a , where 



Ja — 



P 



a 



(85) 



Above ja, the wall velocity X has an oscillating component, while the wall reaches a steady 
motion below j a . The time-averaged velocity is given by 

, sgn[(P-|)i] 



X 



a 



J + 



1 + a 2 



P 



ft 



1 



(86) 
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Fig. 9. (Color online) Threshold current j c plotted as function of pinning frequency f2 = \/2Vo for 
a = 0.01, P = 1, and several values of (3. 




Fig. 10. (Color online) Wall velocity as function of current density for VL = 0.5 and a = 0.01. A jump 
is seen at j = j c . Crossover from (I-a) j c ~ tl/(3 to (I-b) j c ~ tl is seen at j3 ~ § ~ 0.1. 



for j > ja, and 

<*> = f * (87> 

for j < j a . 
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11. Numerical Simulation 

The above results are for a one-dimensional (ID) rigid wall, which would not be the case 
in real experiments (in, e.g., a thick wire of width L± > A). Nevertheless, rigid ID wall 
description seems quite good, if we compare with numerical simulation carried out on realistic 
sample geometries. 22,25 The simulations are based on Landau-Lifshitz-Gilbert equation with 
spin-transfer torque and the /?-term included, without assuming rigid nor ID. It was found 
there firstly that the wall speed is correlated with the appearance of hard-axis component 
of spin (i.e., structure like vortex core) inside the wall. 22 In fact, during slow motion of the 
wall, a vortex core is nucleated, and then the wall is accelerated. In due course, the core is 
annihilated, emitting spin waves, and then the wall slows down and sometimes stops. This 
oscillation of wall speed synchronized with creation and annihilation of vortex core are the 
same as predicted in rigid ID case (Eqs. (71) and (72), where vortex core is simulated by (j)Q 
and periodic modulation of wall speed is represented by sin2(/)o(t) term in the velocity. Thus, 
somewhat surprisingly, deformation and details of spin structure does not affect the dynamics 
in an essential way, resulting in quite a similar wall speed as a function of current (Fig. 10 and 
Fig. 2 of Ref. 22 and Ref. 25 ). This is because of adiabatic wall, where the spin-transfer torque 

3 

tends to flow spin structures at the same velocity j^gPj irrespective of spin structure. 

Effects of deformation would be significant in the presence of extrinsic pinning, and will 
affect the threshold current. 

12. Recent Experiments 

12.1 Metallic systems 

So far experimental results on metallic samples all show threshold currents of order of 
10 12 [A/m 2 ]. If we use K±/ks ~ 0(1) [K] estimated experimentally, 21 ' 57,58 the observed 
threshold is orders of magnitude (10~ 2 — 1CP 1 times) smaller than the intrinsic threshold, 
j l c . For instance, a sample of Yamaguchi 21,57 ' 58 showed j c = 1 x 10 12 [A/m 2 ]. The anisotropy 
energy is estimated to be K±/ks = 2.4[K], and using S ~ ^, a ~ 2.2A and P ~ O(l), 
we obtain j l c = 5.8 x 10 13 [A/m 2 ], i.e., j c /jc ~ 0.02. The observed low threshold currents 
in metals thus should be regarded as due to an extrinsic pinning in regime I-a) or I-b). 
Actually, direct evidence excluding intrinsic pinning in permalloy wires so far was given by 
Yamaguchi et al. 59 They prepared permalloy wires with different geometry, and relaized dif- 
ferent perpendicular anisotropy energies fx-fd — (0.1 — 7.6) x 10 5 J/m 3 , which corresponds 
to K±/ks — 0.03 — 2.4K (per 1 spin). The intrinsic pinning, Eq. (83), predicts then threshold 
current of 5 x 10 11 — 4 x 10 13 A/m 2 . In contrast to this much difference in the predicted values 
of the intrinsic threshold current, experimental values of threshold for these samples do not 
vary so much, (3 — 8) x 10 11 A/m 2 , and are smaller than the predicted intrinsic threshold by 
factor of 2 to 100. Besides, data by Yamaguchi et al. indicate that these experimental values 
do not scale with K±, although there is a weak dependence on K±. Therefore the observed 
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threshold would be of some extrinsic origin. 

Let us first try to explain experimental result 57 assuming regime I-a). Assuming £ ~ A, the 
pinning potential is estimated from the measured depinning field B c = 0.01 — 0.1 [T] as V$ = 
0.34 x (10~ 2 ~ 10- 1 ) [K] = 4.7 x (10~ 26 ~ 10" 25 ) [J], i.e., |- = 1.4x (10~ 3 ~ 10~ 2 ), and so 
j c la ' = (0.21 ~ 0.67) x j l c . This value is still too big to explain the experimental value. Velocity 
jump is estimated as 27 Af Ia ) = — x 839 [m/s], so extremely small (3 (— ~ 4 x 10~ 3 ) is required 
to explain the experimental value of Av ~ 3 [m/s]. 21 If we assume regime I-b), the threshold is 
ic Ib) = 

tot x 2.8 x (10 3 ~ 10 2 ) x j l c . Experimental value could be reproduced if /3 = 0.1 ~ 1. 
But such large value of (3 cannot be explained within the current understanding that (3 arises 
from either non-adiabaticity 9,25 or spin relaxation. 24 Instead, Av cannot be explained by use 
of the above Vo assuming I-b), as it predicts too large value of Av lb ^ = 10 3 [m/s]. Thus, 
honestly, none of the above predictions based on rigid ID wall neglecting temperature rise 
due to heating are successful in explaining experimental result of metals quantitatively. 

It might be crucial to treat the wall as a non-rigid, non-planar object, in particular con- 
sidering the sample width larger than lOOnm. In fact, direct observation of the spin structure 
indicates that the wall is quite deformed upon motion. 60 " 63 It was shown 60 that the initial 
state is not a planar wall but more like a vortex (called a vortex wall), which is the case in 
film or wide wires, and vortex wall moves by applying a current pulse of 2.2 x 10 12 A/m 2 , and 
that the wall is deformed to be a transverse one after some pulses. What was quite interesting 
there is that while vortex wall moves easier, the transverse wall does not move at the same 
current density. Thus the experimentally observed wall motion in wide metallic wires would 
be that of vortex wall, and so simple theory assuming ID rigid wall may not directly apply. 
However, as we discussed, non rigid and non planar nature does not seem essential if we com- 
pare the results to those of numerical simulation. 25 Threshold current of vortex wall obtained 
in simulation is still too large compared with experiments. 64 

There are some possibilities to resolve the disagreement. Most probable one would be the 
heating effect by current. Estimate of Vq by use of experimental B c could be an over estimate 
if effective barrier height Vb is greatly reduced by heating under current, while such heating 
does not occur under static magnetic field. Let us estimate the pinning potential which gives 
the experimental value of j c . Assuming regime I-a), experimental value of j c /f c = 0.02 is 
reproduced if fi = = 1.3 x 10 -5 , which corresponds to V = 3 X 10 _5 [K]= 4.5 x 10~ 5 [T]. 
This is two orders of magnitude smaller than the value extracted from B c . For I-b), we have 
fj, = P x 10~ 2 . From the experiment, Av/(^-j c ) = 3[m/s]/67[m/s]= 0.05. This value is equal, 
for regime I, to ^, so f3 = 5 x 10~ 4 if a = 0.01. So in case I-b), fi = 5 x 10~ 6 . Thus, assuming 
either regime I-a) or I-b), the experimental results could be explained by an extremely weak 
pinning potential, ^5- = 10 -6 ~ 10~ 5 . 

Heating effect in metallic samples has indeed been found to be crucially important. 57 Use 
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of short pulsed current of ns order could be useful in avoiding heating. Sub ns pulse was 
reported to be quite efficient in driving the wall at low current density of ~ 10 10 A/m 2 . 65 This 
could be due to the fact that damping does not affect much for such short timescale. 

Quite recently, Dagras et al. measured the temperature dependence of the threshold cur- 
rent and found that it decreases at low temperatures, for instance, from 2.4 x 10 12 A/m 2 at 
T ~ 170K to 1.9 x 10 12 A/m 2 at T ~ 100K to. 66 Dissipation of spin-transfer torque by spin 
waves was suggested as a possible explanation, but theoretical study is yet to be done. 

12.2 Thin wall 

Quite an interesting result was obtained recently by Feigenson et al. 67 in SrRu03, an 
itinerant ferromagnet with perovskite structure. The current density needed to drive wall was 
5.3 x 10 9 A/m 2 at T = 140K and 5.8 x 10 10 A/m 2 at T = 40K. Small threshold current at 
140K would be due to reduction of magnetization close to Tq = 150K. The threshold current 
is about 2 orders of magnitude smaller than in other metals. This high efficiency would be due 
to a very narrow domain wall, A ~ 3nm, as a result of very strong uniaxial anisotropy energy 
(A') corresponding to a field of 10T. They defined a parameter determining the efficiency 
as a ratio of depinning field and depinning current density, A = B c /j c . Their results were 
A = 1CP 12 Tm 2 /A. They compared this value with threshold current of extrinsic pinning, 9 
given by Eqs. (80) and (81). Using ~ 0.5 x lO" 11 [Tm 2 /A] and S ~ 3/2, we see that 

A = 10 -12 Tm 2 /A is realized if (3 ~ 0.5. This value would be too large if interpreted as 
due to spin relaxation. Using the measured resistivity of domain wall, the non-adiabatic force 
contribution to (5 was estimated and the result of j c was of similar order as observed ones 
but with discrepancy of factor of around 6 at low temperature (Fig. 4(a) of Ref. 67 ). This 
discrepancy seems not so bad considering crude rigid and planar approximation of the wall. 
There is another extrinsic pinning threshold, jc . If we use this expression, A = 10~ 12 Tm 2 /A 
is obtained if K ± ~ n B B c = 4 x (10~ 3 - 10~ 2 )K (per site). 

12.3 Magnetic Semiconductor 

Beautiful experiments were carried out at low current in ferromagnetic semiconductors by 
Yamanouchi et al. 68 ' 69 They fabricated a well structure of 20//m width made of GaMnAs with 
different thickness, which determines the ferromagnetic coupling and transition temperature, 
and trapped a domain wall. The wall position was measured optically after applying a current 
pulse, and the average velocity was estimated. The current necessary was ~ 4 x 10 9 A/m 2 , 
which is 2-3 orders of magnitude smaller than in metallic systems. This is due to the small 
average magnetization, S ~ 0.01, carried by dilute Mn ions, and small hard-axis anisotropy 
Aj_. 69 The obtained velocity was rather high, ~ 22m/s at j = 1.2 x 10 10 A/m 2 . This velocity 
is consistent with the adiabatic spin-transfer mechanism, and the threshold appears to be 
consistent with intrinsic pinning mechanism 9 with anisotropy energy obtained from band cal- 
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culation. However, there are some puzzles. First, the theory of intrinsic pinning 9 and adiabatic 
spin transfer does not take account of strong spin-orbit interaction in semiconductors. So the 
agreement with these thoeries might be a coincidence. 

The second puzzle is the validity of using purely adiabatic theory. In fact, quite a large 
momentum transfer (force) is expected from the wall resistance, R w = 1S7, 70 corresponding 
to P > 1 in terms of /3. 70 

The other puzzle, which was solved just recently, is the temperature dependence of wall 
velocity. The observed velocity scaled as Inv ~ — (Tc — T) 2 j~ 1 / 2 , similar to the creep behavior 
under magnetic field, 71 but this fractional power of j has not been explained in the current- 
driven case. A simple theory of thermal activation assuming rigid wall under the spin-transfer 
torque predicts different behavior, In v ~ j /T, 72 and thus creep motion would be essential in 
the experiment by Yamanouchi et al. Sucessful explanation of creep behavior was just recently 
done by Yamanouchi et al., by taking account of growth of (fr at the pinning center. 

Nguyen et al. studied theoretically the domain wall speed in magnetic semiconductors 
based on the 4-band Kohn-Luttinger Hamiltonian. 38 It was shown there that the wall speed 
can be enhanced by the spin-orbit interaction by a factor of 10 3 — 10 4 due to the increase 
of mistracking, hence reflection, of conduction electrons. This could be useful for efficient 
magnetization switching. 

12.4 Excitation of wall 

Time-resolved study of excitation of wall provides rich information on the wall character 
and driving mechanism. 

Under AC current, domain wall shows another aspect not seen in DC case. AC current 
can drive domain walls quite effectively at low current if the frequency is tuned close to the 
resonance with the pinning frequency. This resonance was realized in recent experiment by 
Saitoh et a/. 55 They applied a small AC current (of amplitude of 10 10 A/m 2 ) in a wire with a 
domain wall in a weak pinning potential controlled by magnetic field. Although the current 
is well below the threshold, the wall can shift slightly as we see below (for about a distance 
of /im, but this would be an overestimate). Under small current, (fro remains small, and the 
equation of motion reduces to that of a "particle" ; 

M W X + — X + M w n 2 X = F{t), (88) 
r 

where M w is the wall mass, r oc a -1 is a damping time, f2 is the (extrinsic) pinning frequency, 
and F(t) is a force due to current. For AC current, I(t) = Ioe lujt , where uj is the frequency, the 
force is given F{t) = ^y-f3 — iPh 2 -^jj , where (3, given by Eq. (73), is from momentum 
transfer and spin relaxation (/3 s f)) arid the last term proportional to uj is from the spin- 
transfer torque. The wall under weak current thus shows a forced oscillation of a particle. By 
measurering the energy dissipation (from complex resistance), a resonance peak would then 
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appear when u is tuned closely around $7. From the resonance spectra, the mass and the 
friction constant were obtained as M w = 6.6 x 10 _23 kg, r = 1.4 x 10 -8 sec. The experimental 
result seems to be well described by the rigid-wall picture, and this would be due to a low 
current density (by factor of 10~ 2 compared to DC experiments on metals), resulting in small 
deformation. What is more, from the resonance line shape, the driving mechanism of the 
domain wall was identified to be the force {(5) rather than the spin-transfer torque. This 
finding was surprising at that time, when adiabatic spin-trasnfer torque was considered as the 
main driving mechanism. The observed force corresponds to the value of j3 ~ 1.5, which is 
too large if (3 arises from spin relaxation /3 s f (/3 s f is considered to be of the same order as a, 
both arising from spin relaxation). If it comes purely from the momentum transfer, the wall 
resistance is estimated to be Row = 3 x 10 [fi], a quite reasonable value. A striking point 
in this experiment is a significant enhancement of the effect of the force due to resonance, 
which made possible the low-current operation. On the other hand, the spin-transfer torque 
is suppressed in the MHz range (as seen from the factor of uj in the spin-transfer torque term 
ofF(t)). 

Quite recently, Thomas et al. 74 succeeded in detecting periodic oscillation of a wall in a 
confining potential by using ns current pulse at j = 6.9 x 10 11 A/m 2 . The motion was consistent 
with the rigid wall description in terms of X and 4>q. Periodic variation of chirality, ^o, of a wall 
was observed in the presence of magnetic field and current pulse of 10ns at 1 x 10 12 A/m 2 . 
The results indicated that the chirality, 0o, plays an important role on the wall propagation, 
as predicted theoretically. 8 ' 9 ' 25 

13. Summary 

We have reviewed theoretical aspects of the current-driven domain wall motion, including 
microscopic derivation of the equation of motion, wall dynamics, and brief discussion on 
experimental results. The effect of current arises from the s-d exchange coupling between the 
local spin and conduction electrons. Treating the non-adiabaticity perturbatively, we have 
derived fully quantum mechanical expression of torques and forces acting on the wall in 
terms of Green's functions. The effect of current on the equation of motion of local spin 
(modified Landau-Lifshitz-Gilbert equation under electric current) was thus obtained. Using 
the results, we derived the equation of motion of the wall. The wall is assumed to be rigid 
and one-dimensional, described by two collective coordinates, position X and angle out of 
easy-plane 4>q. Spin-transfer torque arising from angular momentum conservation was shown 
to contribute to wall velocity, and spin relaxation and non-adiabaticity were shown to work 
as a force on the wall, which induces 4>q. Solving the equation of motion, we found that there 
is a threshold current to drive the wall arising from hard axis magnetic anisotropy energy K± 
and/or extrinsic pinning potential Vq. Threshold current is determined by K± in the intrinsic 
pinning regime, by Vq, K± and force from the current if in extrinsic pinning regime. Our 
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results would be useful in realizing domain wall motion at low current. 

Our formalism can also be applied to describe general spin structures and dynamics under 
current. Extension of our method to first-principle calculation would be useful in realizing fast 
and efficient switching of magnetization by current. 
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